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Abstract 

A model-based LQR method for controlling vibrations in cylindrical shells is presented. 
Surface-mounted piezoceramic patches are employed as actuators which leads to unbounded 
control input operators. Modified Donnell-Mushtari shell equations incorporating strong or 
Kelvin- Voigt- damping are used to model the system. The model is then abstractly formulated 
in terms of sesquilinear forms. This provides a framework amenable for proving model well- 
posedness and convergence of LQR gains using analytic semigroup results combined with 
LQR theory for unbounded input operators. Finally, numerical examples demonstrating the 
effectiveness of the method are presented. 
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1 Introduction 


The use of shell models to describe structural dynamics is pervasive in applications ranging 
from noise reduction in aircraft to flow control in flexible pipes. While general shell equations 
can be used in a variety of geometries, they all share the property that component displace- 
ments are coupled due to the geometry. This leads to significant challenges when developing 
appropriate models and approximation techniques, and constructing effective controllers. 

In this paper, we consider cylindrical shells due to their prevalence in applications. Con- 
trol is provided by piezoceramic patches bonded in pairs to the surface of the shell. These 
transducers provide significant actuating capabilities due to the piezoelectric effect in which 
input voltages generate strains in the patches. Utilization of the converse piezoelectric effect 
(strains produce voltages) also allows the patches to be employed as sensors. When combined 
with their light weight, space efficiency and reasonable cost, these propert ies make the patches 
highly effective control elements in a variety of applications. From a mathematical perspective, 
the use of surface-mounted piezoceramic patches leads to unbounded control input operators. 

Experimental work has already demonstrated the potential for success when employing the 
patches as actuators in applications involving cylindrical shells [8, 13]. However, these initial 
investigations have not. in general, utilized the full potential of the patches due to limitations 
in hardware, models, approximation techniques and control laws. For example, a common 
means of calculating control gains is through the use of modal expansions [10]. However, 
closed form expressions for the modes can be determined only for a limited set of models with 
severely restrictive boundary conditions. The use of incorrect modes when calculating control 
gains can lead to loss of control authority and possible controller instabilities. Hence for most 
ph ysical shell models, modes must first be accurately approximated if modal methods are 
employed for control design. 

In this paper, we present a model-based method for controlling shell vibrations. For sim- 
plicity, the Donnel 1-Mush tari shells equations with Kelvin- Voigt damping are used as a model 
(the assumption of strong or Kelvin-Voigt damping is reasonable and typical for many shell 
materials such as aluminum). The methods are general, however, and can be applied to 
higher-order models (e.g., Byrne-Fliigge-Lur’ye model) if the application warrants. A general 
Galerkin method based on splines is then used to discretize the infinite dimensional system (see 
[7] for details regarding the numerical method and a comparison with finite element methods 
for shells). Through the choice of basis, the method is constructed to be flexible with regard to 
the boundary conditions and material nonhomogeneities which arise in typical applications. 
Furthermore, development of the model and approximation method in terms of a weak or 
energy formulation facilitates consideration of the distributional derivatives which arise when 
including patch contributions in the model. This provides a setting suitable for direct simu- 
lations and control design as well as computation of frequencies and modes for the shell. 

The model and approximate system are then employed in an LQR full state feedback 
theory to obtain feedback gains and, ultimately, controlling voltages to the patches. While full 
state measurements are not available using current instrumentation, and hence the techniques 
cannot directly be implemented in experiments, they provide an important first step in the 
design of effective compensators based on state estimates calculated using a limited number of 
observations (see [5]). The consideration of the LQR performance also illustrates properties of 
the system and model-based control techniques and facilitates investigations regarding issues 



such as patch number and configuration. Finally, the consideration of the problem provides a 
step toward the development of model-based controllers for fully coupled structural acoustic 
and fluid/structure systems involving cylindrical shells. 

The strong and weak forms of the Donnell- Mushtari shells equations are outlined in Sec- 
tion 2. In presenting this model, care is taken to include both passive (material) and active 
(actuator) contributions due to the patches. An abstract form of the model, based on sesquilin- 
ear forms, is also presented. This provides a natural setting to prove model well-posedness and 
convergence properties of the LQR control law. LQR full state feedback laws for systems with 
no exogenous force or forces which are periodic in time are presented in Section 3. In the former 
case, convergence of the approximate suboptimal gains to the optimal gains for the infinite di- 
mensional system is proven using analytic semigroup theory in combination with LQR results 
for unbounded control input operators. A Fourier-Galerkin method for approximating the sys- 
tem dynamics is outlined in Section 1, and the effectiveness of the LQR method for periodic 
forces is demonstrated through a numerical example in Section 5. This example demonstrates 
that through the use of the model-based methodology with general Galerkin approximations, 
significant attenuation in shell vibrations can be obtained using piezoceramic patches. 


2 PDE Model 

The system under consideration consists of a thin cylindrical shell with surface-mounted piezo- 
ceramic patches. It is assumed that the patches are mounted in pairs with edges aligned with 
the circumferential and longitudinal axes of the shell. The edges of the shell are taken to be 
fixed in accordance with common experimental clamping techniques. 

To specify the geometry for the corresponding model, we consider the longitudinal direction 
to be aligned along the .r-axis as depicted in Figure 1. The displacements of the middle' 
surface in the longitudinal, circumferential and transverse directions are denoted by u, e and 
uk respectively while the length, thickness and radius of the shell are denoted by /, /?, /?. The 
region occupied by the middle surface is denoted by Fq* Finally, the shell is assumed to have* 
mass density p. Young s modulus E, Poisson ratio is, Kelvin-Voigt damping coefficient cp and 
air damping coefficient //.. 



x 

Figure 1 . Thin cylindrical shell with surface mounted piezoceramic patches. 
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Actuator and/or sensor capabilities are provided by ,s pairs of surface-mounted piezoce- 
ramic patches. It is assumed that all the patches have thickness h pe , Young’s modulus E pe , 
Poisson ratio i' P , and Kelvin- Voigt damping coefficient ('n,,, ■ Furthermore, it is assumed that 
the glue bonding layer provides negligible contribution to the structural dynamics (the reader 
is referred to [3, 5] for details concerning the incorporation of differing patch characteristics 
and bonding layers in the ensuing models). The region covered by the i th patch pair, with 
edges at xu, x- 2 i, On, 621 , is delineated by the characteristic function 

. m f 1 , x u <x< x 2i , B u <0 < 0 2t 

Xre,(*,0) = < 

I U , otherwise . 


The indicator function S p(t (.r, 0) = Si^{x)Si, 2 ( 0 ) , where 

1 , 6 < (Ou + 0 2 i)/2 

0 , 0 = (8 {l + &2i)/2 

— 1 . 6 > (0u T 0 2 i ) /2 


1 , x <(x u + x 2l )/2 

XiA-r) = \ 0 , x = {x u + x 2 i )/ 2 

— 1 , X > (X \ 1 + X'2 1*)/2 


$ 1 , 2 ( 0 ) ~ { 


delineates the sense of the forces generated by the i th pair. The symmetry of the function 
arises from the property that for homogeneous patches having uniform thickness, equal but 
opposite strains are generated about the point 0 ^ — ((.rj, + a: 2 *)/ 2, ( 0 ]t + 0 2 l )/ 2 ). 


2.1 Strong Form of the Modeling Equation 

We consider here the modified Donnell-Mushtari equations 
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as a model for the thin shell dynamics. As detailed in [3, 5, 9], these equations are obtained 
through force and moment balancing with only low order terms retained. Here M r ,M 9 ,M 9x 
and M xe are internal moments while N x , N$, Ng x and N t6 denote internal force resultants. 
External surface forces are denoted bv q r , qe, q n whereas the external resultants (line moment s 
and forces) generated by the i th patch pair are designated by ( M r ) pfl , ( M e N r ) pc , , ( No )^., . 

Expressions for the internal force and moment resultants are derived under the assumption 
that stress is proportional to a linear combination of strain and strain rate. This yields a. 
model which incorporates Kelvin- Voigt or strong internal damping. As detailed in [3, 5], the 



resultants N x ,N x e, N $ x , M x , M x $* Mq x derived under this assumption are 
Nx = -j— ■ 2 + 1/60 ) + 1 JVt r ) \ pe. ( X 1 0 ) 
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where the constant «;* = (h / 2 + h p < )'"* — // 3 /8 results from integration through the thickness of 
the patch. Expressions for the resultants N$ and A1$ can be obtained by replacing s x .6q, k x , kq 
in the expressions for N x and M x by So, t x , k$, k x . respectively. The midsurface strains and 
changes in curvature for the Donnell-Mushtari model are 
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Note that for the undamped shell which is devoid of patches, the resultant equations (2.2) 
reduce to the classical Donnell-Mushtari expressions 
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To characterize the external or active patch contributions, it is typical to start with the 
assumption that the strains generated by a patch are proportional to the applied voltage 
[3]. Since differing voltages can be applied to the outer and inner patches in the pair, we 
will differentiate between the two with V tl (t) and \h{t) used to denote the voltages to the 
outer and inner patches in the i tb pair, respectively. The proportionality constant relating 
the generated strain to the input voltage is designated by d 31 . As detailed in [3], the total 
external moments and forces generated by the patches are 
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where u> = (h/2 + h pf ) 2 — /i 2 /4 and u 3 = (hj 2 + h pi ) ! — /i 3 /8. When substituted into (2.1), 
the expressions (2.4) provide the input from the patches when voltages are applied. 

Finally, the fixed-edge boundary conditions 


die 

u = v — w = = 0 

ox 


x = 0, f. 


(2.5) 


are used to model the end behavior of the shell. These boundary conditions are appropriate 
for experimental setups in which heavy endcaps prevent edge movement. Note that alternative 
boundary conditions such as simply supported or “almost fixed ” (se e [4]) can be employed if 
edge movement is suspected. 


2.2 Weak Form of Modeling Equations 

The strong form (2.1) of the modeling equations requires first and second derivatives of the 
moment and force resultants. As noted in (2.2) and (2.4), both the internal and external 
moment and force resultants are discontinuous due to the piezoceramic patches. Hence formal 
analysis and approximation using the strong form of the modeling equations lead to difficulties 
due to differentiation of Dirac distributions. 

To alleviate these difficulties, it is advantageous to consider a weak form of the modeling 
equations which can be derived from Hamilton’s principle (energy considerations). While 
equivalent to the strong form under suitable smoothness assumptions, the weak form provides 
a. more natural setting for analysis and approximation. 

The state variables for the problem in second-order form are taken to be y = (u, i\ w) in 
the state space H = I 2 { ID) x /. ‘ ( l’y ) x I 2 ( 1 0 ) . For the fixed-edge boundary conditions (2.5), 
the space of test functions is taken to be V = / ( I t , ) x 7/ ( J ( T 0 ) x //^(F 0 ) where 

//d(F„) = {»/ € 7/ 1 ( T 0 ) | ?/(0, 9 ) = = 0} 

tfo(Fo) = in € // 2 (r o )|j/( 0 ,tf) = ,J,(O,0) = v((J)) = yA<-0) = 0 } . 



For <!> = (u,v,w) and 'k -- (//i . ?/ 2 - ) -* the H and l inner products are taken to be 

(<t>. = / phuTpd-, + / phvT) 2 d-) + / phwT^dp/ 

./ 1 o J To ^ 1 o 

and 

/•:/» 


<(£\ £■„,)*, *> v . = 


dip 1 , (hit 

If ' + " £,) a7 + 2fl (1 “ ^ w 


^ 2 hpphpf 


i = i 4 *V 




, 1 , v ^77* 

( ' r + Up, ' 9) lh + w {l ~ Vpr)Sre Iff 


+ 


Eh 

do I 


, . _ , dip 1 (h/2 

l ^ + UZi} ^ + w ([ ~' ,] - x6 lh 




,=> ■ - V 


, . , d m l () i/ 2 

(■" + pt:j ) m + -2/(“ ( ' v)c ^ ox 


+ 


Eh 

./[„ (I - '-' 2 

h 2 


i , V , <) 2 m 

-(eg + - J^(Kr + VKg ) “Y 


, , /P'/3 t> 2 j , , - 


^: >(| 2 -"^) V, "' (I ’ >) 
i=l v P'' 


3/* 


/? 


( - 0 T + Vp f Kq ) ^^2 


«:i , , r/ i M v ^ 2/ /3 

r2 \ 0t~ p, jj ml R \ P ej d ^ d0 


dr 


where €&* 5^, k j: . hq, t are defined in (2.3) and r/y = UdOdr. The dependence of the inner 
product on tin 1 Young’s moduli is explicitly included in tlx 1 definition to provide a notation 
for defining analogous damping expressions later in this work. 

The weak form of (2.1), as derived in [5] from energy principles, is given by 
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for all = (r/i,?/ 2 ,r/ 3 ) G V. A comparison between (2.6) and (2.1) illustrates that in the 
weak form, derivatives are transferred from the discontinuous resultants onto suitably smooth 
test functions. This alleviates the difficulties associated with the discontinuities and reduces 
smoothness requirements on approximate solutions. 


2.3 Abstract Formulation 


To define appropriate sesquilinear forms, we group stiffness components separately from damp- 
ing components. To this end, we define a, : V x V' ->€,/= 1,2 by 


*.(*,#) = {(E, E pe )*,*) v 

,t 2( < 1 ) . 'b ) = ((e/p,CD p J4», '1'^. + ' 

Note that ((c/y, e/y ; „ )<1>, 'P)v differs from ((E.E pc )< h, V l , ) v - only in that Young’s moduli are 
replaced by Kelvin- Voigt damping coefficients. It can be directly verified that the stiffness 
form <T| satisfies 
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(Symmetric) 
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Remark 1. The symmetry of ey is dependent upon the choice of shell model and ultimately 
reflects the Maxwell-Betti reciprocity theorem. W hile the Donnell-Mushtari model yields a 
symmetric sesquilinear form < 7 ,, other models such as the Timoshenko shell model will not 
yield a symmetric form. 


To represent control contributions, let V = III s denote the Hilbert space of control inputs 
and define B G C(U,V*) by 
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for f 6 l . Here {•.•) v , r denotes the usual duality product. Finally, with the definition 
<1 = ( • / pl > , qe. cy w ], we can write the weak form (2.6) in the abstract variational form 


(y(t)-*)v.v + *20/(1),*) + <Ti(y(t),V) = (Bu(t) + f j(t),y V) v . v . (2.8) 


< 



To pose the problem in a first-order form amenable for control applications, we define the 
product spaces H = V x // and V = V x V with the norms 

I ( 0 1 , 02 ) | = |0i lr + |02lw 

|(0ii0'i)lv — 1 0i lv + | 02 |r • 


The state is taken to be z(t) — (y(t).y(t)) G H. Finally, the product space forcing terms are 
formulated as 


</(/) = 
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Bu(f) — 
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Bu(t) 


(2.9) 


The weak form (2.8) can then Ik* rewritten as 


(~(/), A) v , v + it(~(/). A) —{/?«(/) + </(/), A) v . v for A G V ^ ^ 

~( 0 ) = = (yo,yi) 

where rr : V x V — >■ d' is given by 

rr(0. r) = - (02, 0i)v- + o-i (01, 02) + ^ 2 ( 02 , k’ 2 ) 

for 0 = (0i,02), 0 = (0i - 02 ) € V. As proven in [b, page 109], rr is V continuous and for 
A > 0, <t(-, •) + A (•. -) v is V-elliptic. From the coiitimiity of a. il follows that one can define 
an operator A G C(V,V*) by <j(T,A) = (*4T, A)y*,v * 

To obtain a strong form of the first-order system which is appropriate for control purposes, 
consider the system operator 


dom*4 — {(0i,02) G G V, A\ 0 \ -|- A 2 O 2 G //} 


A = 


0 / 

-A t -A 2 


(2.11) 


with Tj , A 2 G £(V, l *) defined by 


(Ai<l>t,<h)v\v ~ fT <( 0 i: 02 ) , i = F2 . 

It should be notated that A is the negative of the restriction to dom.4 of A € £(V.V*) so 
that <x(T, A) = (-.4T. A) w for T € domAA € V. 

A strong form of the abstract system model is then given by 

~(t) = Az{t) + Bu{t) + g(t) ill V* = V x v m 

( 2 . 12 ) 

=(0) = c o . 

'The rigorous equivalence of solutions is established through the following theorems. 


Theorem 1. Under Hypotheses (U1)-(H5) on rr x and rr-i. A generates an analytic semigroup 
T(0 on V,H and V*. In terms of this semigroup, the representation 

-(/) = T (/)- 0 T f T(i -s)[Bu(s) + A,)]ds ( 2 . 13 ) 

■to 

defines a mild solution to (2.12) for c () G V* and Bu + y G (0, 7 ); V*). Furthermore, this 
semigroup is (uniformly) exponentially stable on V. H and V*. 
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Theorem 2. Let z sg denote the semigroup solution to (2.12) given by (2.13) and let r„ nr denote 
the weak solution to (2.8). Under hypotheses (HI )-(II5), it follows that, c S5 (c 0 , T) = z var (z 0 ,f) 
for z 0 £ "H and T = Bn + g € L 2 ((0, T); V*). 

Following the convention of [14], we will use the same notation for the semigroups defined 
on Vh'H and V* since each semigroup is an extension or restriction of the others. Note that 
dom«4 defined in (2.11) is actually dom^.4, the domain of A as a generator of T(t) in H. As 
detailed in Lemma 3.6.1 and Theorem 3.6.1 of [14] (see also Section IV. 6 of [12] and Chapter 2, 
Theorem 5.2 of [11]), the property that A generates an analytic semigroup on V, H and V* 
results from the continuity and V-ellipticity of a. The exponential stability of T(t) on H 
for second-order systems with strong damping is demonstrated in [1] while the exponential 
stability of T{t) on V and V* in this setting is proven in Lemma 3.3 of [2]. Finally, Theorem 2 
is a reformulation of Theorem 4.14 of [5] and details can be found therein. 

3 LQR Control Problem 

In the last section, the PDE system modeling the dynamics of the thin shell with surface- 
mounted piezoceramic actuators was written in the abstract first-order form 

z{t) = Mt) + Bu(t) + g(l) 

~(0) - ~o 

in V*. In this section, LQR control results for both the original infinite dimensional problem 
and approximating finite dimensional problems will be discussed. Two cases will be consid- 
ered, namely when g = 0 and g is periodic in time. In both cases, it is assumed that state 
observations in an observation space Y have the form 

z oh (t) = Cz(t) (3.1) 

where C £ C(H , V) is bounded. The assumption that C- is bounded is made to simplify the 
exposition and the reader is referred to [2] for arguments pertaining to the case of unbounded 
observation operators. 

3.1 No Exogenous Input 

For the case in which g = 0. the infinite horizon problem concerns the determination of a 
control u which minimizes the quadratic cost functional 

J(u,z 0 ) = jH {\Cz(t)\i- + \K l 'Mt)\l}dl (3.2) 

subject to 

i(l) = Az(t) + Bu(t) 

5(0)= C 0 . 

I he positive, self-adjoint operator 7? is used to weight various components of the control. 
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As detailed in [2, 5], if (A,B) is stabilizable and (A,C) is detectable, then the optimal 
control minimizing (3.2) is given by 

u(<) = 

where II solves the algebraic Riccati equation 

{A * n + r ia - n bk- 1 b*ii + cc)z = o for ail ~ g v 

and z(t) ~ 5(0-o. Here 5(/) is the closed loop semigroup generated by *4 - 57? -1 5*n. 

For implementation purposes, it is necessary to define an approximate system and controls, 
arid determine convergence criteria for these approximate controls when fed back into the 
infinite dimensional system. The approximations are considered in a Galerkin framework 
with trajectories evolving in the finite dimensional subspaces V A C V C H. It is assumed 
that the approximation method satisfies the standard convergence conditions 

(H1N) For any c £ V, then' exists a sequence c A £ V A such that |~ — r A |v — > 0 as N — > oc. 

The finite dimensional operators and approximating system an 1 t hen determined as follows. 
The operator *4 A : V A — > V A which approximates A is defined bv restricting a to V A x V A ; 
this yields 

(-A* T. A) = cr(T, A) for all T. A G V N . (3.3) 

For eacli A', the fb semigroup on V A which is generated by *4 A is denoted bv T A (/). The 

control operator is approximated by 5 A G £{ ( . V A ) given bv 

A^ — (u„B* A)y for all u G P , A G V A (3.4) 

while CA denotes the restriction of the observation operator C to V A . Finally, we let / ;A 

denote the usual orthogonal projection of H onto V A which by definition satisfies 

(i) p v y e v jV for T £ H 

(ii) (p v T - T, A) = 0 for all A G V A . 

This projection can be extended to G £( V*, V A ) by replacing t he 'H-inner product (T, A ) y 
by t he duality product (T, A) v . v and considering T £ V*. 

The approximate problem corresponding to (2.10) with <j = 0 can then be formulated as 

^(r v (/). \) w • rr(, v (/).A) = A) w for all A € V v 

: V (0) = P N z 0 . 

This has t he solution 

c v ( / ) = T v ( / ) P N Co + [' T N (t-#) P* B s u ( .s )ds . 

Jo 
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The following theorems taken from [2, 5] can he used to establish the convergence of the 
approximate gains to their infinite dimensional counterparts for certain classes of shell models 
(see specifically Theorem 7.10 and Lemma 7.13 of [ 5 ]). 

Theorem 3. Assume that the injection i : V H is compact. Moreover, suppose that the 
damping sesquilinear form can be decomposed as a 2 — &<J\ + dy, for some S > 0 , where the 
continuous sesquilinear form a 2 satisfies for some A G IR 

Re d 2 (d>, <f>) > y ~ A|0|/y for all <j> G \ . 

Finally, suppose that the operator A"[ l A 2 , where \ 2 G £( V . V *) is defined by (^A 2 (f>, = 

d 2 ( 0 , //). is compact on V . 

If for some u: G Hi and M > 1 . T(t) satisfies 

I T(t)\c(n) < A/c u;/ , / > 0 , 

then for any e > 0 there exists an integer N* such that for N > Ay, 

|T Ar (/)P N | £( ^) < , / > 0 

for some constant M > 0 independent of N . 

Theorem 4. Assume that the injection i : V <—» // is compact. Let the sesquilinear form o 
associated with the first-order system (2.10) be continuous and V-elliptic. Assume that the 
operators AXX of ( 2 . 11 ), (2.9), (3.1), respectively, satisfy: (AX) is stabilizable and (AX) 
is detectable where B G C(IKV*) is unbounded and C G £[%,}') is bounded. Consider an 
ap])roximation method which satisfies (H1N). Finally, sup])ose that for fixed A'o and N > A' a , 
the pair (*4 A ,/? A ) is uniformly stabilizable and (*4 A ,C A ) is uniformly detectable. 

Then for N sufficiently large, there exists a unique non negative self-adjoint solution 0 A G 
£(V*, V) to t he N th approximate algebraic Riccati equation 

A N *n N + n N A N - II N B N K- l B N 'll N + C N *C N = 0 

in V A . There also exist constants A /3 > 1 and u! :i > 0 independent of N such that <S A ( / ) = 
t iA N -&xiz- l B Nm n N )t satisfies 

S N (t)\ vN <M; i( -^ , / X). 

Moreover, the convergence of the Riccati and control operators 

n A P '' ~ A Hr in V for every : t F 
B N’Y[NpN _#*n| -> 0 , 

as A r — >oc. is obtained. 
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Example 1. We consider in this example a shell with constant parameters p, E, v, cp. Such a 
case would arise if modeling a homogeneous shell or a shell in which the variance of material 
properties across regions with actuators is negligible. The sesquilinear forms for this model 
are specified in (2.7). Due to the constant coefficients, cr 2 can be written as a 2 — Scr { + a 2 
where 5 — ^ and T) — p Jj- wrfad^. It follows immediately that 

H((Ti{0,4>) = f.1 / 0 2 d~j > — — |<2>| V- 

for all cj) G V . The boundedness of the operator A 2 generated by a 2 follows directly from 
the boundedness of a 2 . Furthermore, it is noted that /If 1 G £(V*,Y ) can be written as an 
operator on V V by If 1 = .If 1 /*/ where the injections i : V c — > //,/'* : // V * are 
compact. Thus A f 1 is comj>act on V which implies that .-If 1 . 4 2 is compact on V since it is 
formed from the product, of compact and bounded linear operators. Finally, the exponential 
stability of T(t ), the stabilizability of (A, B) and tlie detectability of (A, B) are guaranteed by 
Theorem 1. The hypotheses of Theorem 3 are then satisfied for this system and one obtains 
uniform bounds on the approximating semigroups. Thu convergence of the Riccati and control 
operators is then obtained from Theorem 4. 

3.2 Periodic Exogenous Input 

A reasonable assumption in many mechanical systems is t hat g is periodic in time with period 
r. t he system to be controlled in this case is 

i(l) = Az{t) + Bu(t) + g( t ) 

(T- r >) 

~(0) = *(r) 


and an appropriate quadratic functional to be minimized is 

•C(«)=f £ {\C~(t)\y + |^ 1/2 «(0|( } ■ 

l Jo 

Note that the periodic exogenous term g can bo used to model inputs such as noise gener- 
ated bv rotating engine' components (e.g., propellers or turbines) or periodic electromagnetic 
disturbances. 

To guarantee the existence of a unique Riccati solution and control for the system (3.5), 
it is assumed that (A,B) is stahilizable and (A.C) is detectable. Furthermore, it is assumed 
that g G L 2 (0,r:H) and that B is bounded. Under these conditions, it is verified in [Oi] that 
ill*' Riccati equation 

.4*11 + UA + IWR-'B'U + C‘C- = 0 

lias a unique solution. Furthermore, if r denotes the r-|>eriodic solution of the adjoint or 
tracking equation 

/■(/) = -[A - BR- l B'\l]*r(t) + II (j{t) 
r(0) = r(r ) 
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and z is the closed loop solution of 


~{t) = [A - BR- l B*n]z(t) - BR~ l B'r(t) + <,(t) 

A 0) = ~(r), 

then the optimal control is given by 

u(t) = --R-'B*[n~(n-r(t.)}. (3.6) 

The LQH theory for this case is less complete than that for systems with no exogenous 
input and is currently limited to bounded control inputs B. The synthesis of the theory for 
unbounded input operators and periodic exogenous forces is currently under investigation. 
The effectiveness of the method is illustrated in the final example of this work. 


4 Approximation Method 

A Galorkin method was used to approximate the solutions n, v, ir to the system (2.6), or 
equivalently, (2.10). The approximating subspaces were taken of the form V A = span{£? ttfc }x 
span{£h A } x spanj/3,^} where B Uk , B Vk , B Wk denote bases for the ■*/., v and w displacements, 
respectively. To exploit the tensor nature of the shell domain F 0 and periodicity in 0 , the bases 
were constructed with Fourier components in 0 and cubic splines in x (see [7] for details). The 
approximate' displacements were then given by the expansions 

A r u 

(t,0,x) = u k(t )B Uk {0, x ) 

k=i 
Mv 

v N (t,0,x) = v k {t)B Vk {Q,x) 

k=\ 

A r v> 

N (tA*) = £m(t)B Wk (0,x). 

k= i 


w 


To obtain a finite dimensional system with matrices corresponding to the finite dimensional 
operators in (3.3) and (3.4), the sesquilinear forms <7\ and <r 2 were restricted to V A . This yields 
the matrix system 


' 1<E 

0 ■ 

■ tf A >) ■ 

■ 

0 

\/ A ' . 


- 

' Af 

0 ' 

’ tf A '(0) ' 

r 

0 

\/ A ' . 

t) A '(0) 

L 


0 

-Af 


y\ 

..A' 






+ 


0 

/i A ' 


[«(/)] + 


0 


<j A (0 J 


where i) A (i) = [■«, (# ), • • • , ujy u , ( t ) . • • • , • • • . tPj*,'..] 7 ' contains the Af — A r u +Af l . + A r „. 

generalized Fourier coefficients. The .s patch inputs are contained in u(t) — [»](/), • • • ,u s (t)] 1 . 
The reader is referred to [7] for details concerning the construction of the mass, stiffness and 
damping matrices A/ A . A a , A a , the inputs /f A . </' V (/) and the init ial conditions y A ' 2/ 2 • 



Multiplication by the inverted mass matrix yields the Cauchy equation 

i N (t) = A N z N {t) + B N u(t) + g N (t) 

=■» = 

where z^ € /? 2A = [i'F v (t)zd* (/)] T * This system forms the constraint equations used in the 
finite dimensional LQR theory discussed in Section 3. 

5 Numerical Example 

We consider here an exogenous force g which is periodic in time with period r — IOGOtt 
( 500 Hz). The distribution of the force was taken to be binomial in the transverse and 
longitudinal directions and centered at {x,0) = (f/2,0) and (xj)) = (f/2,7r) as depicted in 
Figure 2. The magnitude of the transverse component q n was one hundred times that of the 
longitudinal component q x so as to model an input consisting primarily of acoustic sources 
located adjacent to (f / 2,0) and (//2, 7r). 

Six pairs of piezoceramic patches of length lem and radial measure 7t/3 were employed 
as actuators. The locations and material properties of the patches along with the dimensions 
and physical parameters for the shell are summarized in Table 1. 

To accommodate the periodic exogenous force </, control inputs to the twelve patches wen* 
computed using the feedback law (3.6). Note that in this formulation, independent voltages 
an* determined for the individual patches. This provides the capability of generating both 
inplane forces and bending moments in the regions covered by tin* patches so that longitudinal, 
circumferential and transverse vibrations can be controlled. 

Time histories of the uncontrolled and controlled shell displacements at the point /q = 
(•]( /4, 7r/32), depicted in Figure* 2, are plotted in Figure 3. The open loop trajectories exhibit 
both a transient response settling into steady state and a beat phenomenon due to the close 
proximity of the driving frequency and natural frequencies for t he shell. At t his observation 




0 - Distribution 
of Normal Force 


Figure 2. Distribution of normal forcing function at 0 — 0 and 0 = tt. Observation lines 
L\ — {(.r,fl)|0 < x < fj) — 7t/ 6}. A 2 — {(;r, 0)|.r = 3f/1.0 < 0 < 2tt } and observation point 
P\ = (3C/4. 7r/32). 


11 



point, all three displacement components are reduced by more than 90% when controlling 
voltages are fed back to the patches. 

To illustrate the spatial attenuation due to the feedback of voltages to the patches, root 
mean square (mis) plots of the uncontrolled and controlled trajectories along the axial line 
L\ and circumferential line L 2 (see Figure 2) are plotted in Figure 4 and 5, respectively. 
For the open loop case, these plots illustrate a standing wave in all three components of the 
displacement (the slight asymmetry in the axial plots is due to the longitudinal input q r ) . 
The figures also demonstrate significant reductions in all three displacement levels, even in 
regions not covered bv patches. This further illustrates the effectiveness through which the 
model-based control law can be used to attenuate shell vibrations. 






2 


3 


re (mis) displacements along 
out rolled). 


]() 








Dimensions 

Parameters 


h — .00127 m 

p = 2700 kg/rn 3 


R = .4 m 

E = 7 . 1 x 10'° N/m 2 

Shell 

( = 1 m 

cj) = 1.17 x 10 r> Nws 

V = Ml 

fi = 58.97 Ns/m 2 


h pe = .0001778 m 

Ppe = 7000 kg/ m 3 
E p( = 0.4 x 10 10 N/m 2 

Patches 

Centers (x,0): (.25,0), (.5, 0), (.75, 0) 

vT 

II 

X 

o 

ca 

> 

> 


(.25, 7T ), (.5. 7r), (.75, 7 r) 

i'pt = .51 


Dimensions: x : 0.1 an, 0 : n/'i 

d-M = 100 x 10 -12 m/V 


Table 1 . D imensions and physical parameters for the shell and patches. 


6 Conclusions 

A model-based LQR method for controlling shell vibrations has been presented here. While 
developed in the context of a modified Donnell-Mushtari cylindrical shell model, the method 
is quite general and can be directly extended to other models and geometries. Under the 
assumption of strong or Kelvin- Voigt damping (a reasonable and typical assumption for many 
shell materials), model well-posedness arid convergence of control gains is obtained using 
analytic semigroup theory combined with LQR results for unbounded input operators. 

The Galerkin method used to approximate the system dynamics utilizes bases constructed 
from tensored Fourier polynomials and modified cubic splines. As discussed in [7], care must 
be taken when developing methods for approximating shell dynamics so as to avoid shear or 
membrane locking. One manifestation of locking is the existence of model dynamics which 
are incorrect ly approximated by the numerical method. The use of a numerical method which 
exhibits locking can lead to a loss of control authority and potential controller destabilization 
if the approximations are sufficiently inaccurate. Further details regarding issues concerning 
the approximation of shell dynamics and convergence properties of the numerical method can 
be found in [7]. 

The numerical example demonstrates the effectiveness of the model-based control method 
for attenuating all three components of the shell displacement in the presence of both transient 
and steady state dynamics. Furthermore, by modeling the global shell dynamics and patch 
interactions through coupled PDF and constructing the control law in terms of these PDE, 
significant reductions in displacement levels throughout the shell are obtained, even in regions 
devoid of patches. Numerical implementation of the LQR method in this manner provides a 
first step toward the development of model-based state estimators and compensators which 
can be experimentally implemented in shell applications. 
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